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Physics  of  Atomic  Frequency  Standards  Based 
on  Stimulated  Resonance  Raman  Scattering 


1.  BACKGROUND 

In  this  report,  we  analyze  an  attractive  possibility  for  achieving  high  precision 
frequency  standards  based  on  stimulated  resonance  Raman  scattering  in  atomic 
and  molecular  beams.  We  have  proposed  the  use  of  this  effect  in  frequency  stand¬ 
ards,  as  well  as  a  first  embodiment  in  a  sodium  atomic  beam  as  an  outgrowth  of 

1  2 

an  ongoing  RADC  program.  ’  The  presentation  will  be  tutorial  in  nature  and  will 
emphasize  the  physics  of  the  interaction  relevant  to  device  operation,  while  avoid¬ 
ing  many  important  engineering  considerations. 

The  standards  being  considered  rely  neither  on  inhomogeneous  magnets  for 
polarizing  or  analyzing  the  atomic  or  molecular  beam  nor  on  hot-wire  detectors  for 
atomic  detection.  Rather  they  employ  lasers  for  state  preparation  and  detection  of 
transitions,  as  well  as  for  the  basic  interaction.  Through  a  stimulated  Raman 
process,  the  frequency  difference  of  two  lasers  striking  the  atomic  or  molecular 
beam  simultaneously  in  the  interaction  region  is  stabilized  to  a  long-lived  atomic 
or  molecular  transition.  The  difference  of  the  two  laser  frequencies,  and  hence 

(Received  for  publication  19  October  1981) 

1.  Ezekiel,  S. ,  Leiby,  Jr.,  C.C.,  Picard,  R.H.,  Willis,  C.R.,  and  Hackel,  R.P. 

(1980)  Laser  Stimulated  Raman  Molecular  Beam  Time  and  Frequency  Stand¬ 
ard,  allowed  U,  S,  patent  application  No.  150,  523,  filed  16  May  i960. 

2.  Thomas,  J.  E. ,  Ezekiel,  S.,  Leiby,  Jr.,  C.C.,  Picard,  R.H.,  and  Willis, 

C.R.  (1981)  Opt.  Lett.  6:298. 
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the  frequency  output  by  the  atomic/molecular  portion  of  the  standard,  may  vary 

from  the  microwave  to  the  infrared  region  of  the  spectrum,  depending  on  spacing 

of  energy  levels  involved  in  the  transition. 

We  concentrate  in  this  report  on  factors  affecting  the  precision  of  the  frequency 

standards.  For  definiteness,  we  envision  stimulated  Raman  scattering  between 

o 

the  hyperfine  levels  (F=l  and  F=2)  in  the  ground  electronic  state  (3  Sp^j)  the 

sodium  atom,  separated  by  the  microwave  (L-band)  frequency  1.772  GHz,  For  the 

moment,  we  ignore  any  subtleties  introduced  by  the  magnetic  substates  (Mp  =  0, 

±1  for  the  F=1  level  and  Mp  =  0,  ±1,  ±2  for  the  F=2  level).  This  transition  is 

the  sodium  transition  analogous  to  the  hyperfine  transition  (F=4  —  F=3)  at 

o 

9. 193  GHz  induced  in  the  6  Sp^  electronic  ground  state  of  cesium  in  the  usual 

cesium  atomic  clock.  To  enhance  the  Raman  cross  section,  we  work  near  optical 

2  2 

resonance  of  the  two  laser  beams  with  the  3  or  ®  ^3/2  exc^tec^  state  of  sodium. 

The  single-photon  emission  between  the  32P  levels  and  the  ground  state  gives  rise 
to  the  very  intense  characteristic  yellow  doublet  D  lines  (Dp  and  Dg)  of  sodium. 

The  history  of  the  Raman  effect  began  fifty  years  ago.  The  existence  of 
scattered  light  with  frequency  shifted  by  an  amount  corresponding  to  the  energy 

difference  between  two  quantum  states  was  discovered  in  1928  by  Raman  and 

3  4 

Krishnan  in  liquid  solutions  and  by  Landsberg  and  Mandelstam  in  solids.  This 

discovery  led  to  ordinary  or  spontaneous  Raman  scattering  spectroscopy,  which 
today  is  used  mostly  as  a  diagnostic  tool  for  studying  the  vibrational  energy  levels 
of  molecules  and  of  optical -branch  lattice  vibrations  in  crystals  and  for  identify¬ 
ing  unknown  samples  or  trace  constituents.  A  spectral  analysis  of  the  scattered 
radiation  reveals  the  existence  of  frequencies  shifted  down  (Stokes  scattering)  or 
up  (anti-Stokes  scattering)  by  increments  equal  to  vibrational  frequencies  of  the 
material  irradiated.  A  second  weak  beam  of  light  applied  to  the  sample  at  the 
Stokes -shifted  (or  anti-Stokes -shifted)  frequency  will  experience  gain.  This  is  the 
phenomenon  of  stimulated  Raman  scattering  first  observed  after  the  invention  of 

5 

the  laser  by  Woodbury  and  others  In  1962,  In  conjunction  with  vibrational  transi¬ 
tions  in  organic  laser  Q-switching  liquids.  Of  course,  it  is  not  necessary  that  the 
initial  and  final  states  of  the  scattering  sample  be  different  vibrational  levels,  as 
they  were  historically.  Any  two  quantum  states  of  the  sample  having  the  same 
parity  will  do. 

In  the  quantum  theory  of  radiation,  Raman  scattering  appears  as  a  two-photon 
process  involving  the.J  •  A  term  of  the  interaction  Hamiltonian,  where is  the 


3.  Raman,  C.V.,  and  Krishnan,  K.  S.  (1928)  Nature  (London)  121; 501. 

4.  Landsberg,  G.,  and  Mandelstam,  L.  (1929)  Naturwtssenschaften  16:557,772. 

51  Woodbury,  E.J.,  and  Ng,  W.  K.  (1962)  Proc.  IRE  50:2367:  Eckhardt,  G., 
Hellwarth,  R.W. ,  McClung,  F.J.,  Schwarz,  S.E. ,  Weiner,  D. ,  and 
Woodbury,  E.J.  (1962)  Phys.  Rev.  Lett.  9:455. 
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polarization  current  density.  Since  this  Interaction  is  linear  in  the  field  vector 
potential  A,  it  causes  transitions  Involving  two  photons  through  the  agency  of  inter¬ 
mediate  states  of  the  scattering  sample  and  the  light  field  which  differ  from  initial 
and  final  states  by  having  only  one  photon  absorbed  or  emitted.  The  Raman  cross 
section  will  involve  a  sum  over  all  possible  intermediate  states.  Energy  need  not 
be  conserved  in  the  intermediate  state,  but  only  in  the  over-all  transition  from 
initial  to  final  state.  However,  the  proximity  of  the  intermediate-state  energy  to 
an  actual  atomic  resonance  leads  to  a  large  enhancement  of  the  Raman  cross  sec¬ 
tion  similar  to  the  enhancement  of  Rayleigh  scattering  near  an  atomic  resonance 
(resonance  fluorescence).  This  enhancement  of  the  Raman  effect  is  known  as 
resonance  Raman  scattering.  In  the  resonance  Raman  effect,  there  is  only  one 
intermediate  state  of  any  consequence  —  the  resonant  one  —  and  the  sum  over 
intermediate  states  reduces  to  a  single  term.  Hence,  one  need  deal  only  with 
three  energy  eigenstates,  the  initial,  the  final,  and  the  one  intermediate  state  of 
consequence.  It  follows  that  the  physics  of  resonance  Raman  scattering  is  inex¬ 
tricably  tied  up  with  the  physics  of  material  systems  having  three  quantum  states, 

g 

or  three-level  systems.  In  Section  2,  we  will  analyze  stimulated  resonance 
Raman  scattering  in  terms  of  the  semiclassical  radiation  theory  of  an  optical  field 
interacting  with  three-level  systems,  referred  to  hereafter  as  three-level  "atoms." 

Before  proceeding,  we  should  point  out  an  important  characteristic  of  reso¬ 
nance  Raman  scattering,  namely,  it  cannot  be  described  as  absorption  of  the  inci¬ 
dent  photon  followed  by  emission  of  the  scattered  photon.  This  was  first  demon¬ 
strated  convincingly,  in  conjunction  with  resonance  fluorescence,  about  1930  by 
7  8 

Placzek  and  by  Weisskopf.  Rather,  one  must  view  the  transfer  of  the  atom 
from  initial  to  final  state,  accompanied  by  destruction  of  an  incident  photon  and 
creation  of  a  scattered  photon,  as  a  single  quantum  process.  In  other  words, 
transfer  of  atoms  from  the  initial  to  the  final  state  is  not  a  two-step  process, 
comprised  of  successive  single-photon  interactions,  but  a  single-step  coherent 
two-photon  interaction.  An  Important  practical  consequence  of  this  is  that  the 
atom  resides  in  the  intermediate  state  for  a  negligible  amount  of  time  and,  con¬ 
sequently,  the  linewidth  of  the  upper  level,  which  can  be  very  large,  does  not 

Q 

adversely  affect  the  width  of  the  Haman  resonance. 

2 

In  sodium,  for  example,  the  3  P  Intermediate  state  is  separated  by  optical 
o 

photon  energies  from  the  3  S  initial  and  final  states,  and  the  3P-3S  transition  is 


6.  See,  for  example,  Beterov,  I.  M..  and  Chebotaev,  V.  P.  (1974)  Progress  in 

Quantum  Electronics.  Vol.  3,  Part  1. 

7.  Placzek,  G.  (1929)  Z.  Phys.  38:585. 

8.  Weisskopf,  V.  (1931)  Ann,  Phys.  (Leipzig)  9:23. 

9.  Schlossberg,  H.R.,  and  Javan,  A.  (1966)  Phys.  Rev,  Lett.  17:1242. 
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strongly  allowed;  consequently,  the  3P  state  has  a  large  natural  linewldth  of 
10  MHz.  On  the  other  hand,  the  initial  and  final  states  are  either  the  atomic 
ground  state  or  levels  lying  close  to  the  ground  state  whose  transitions  to  the 
ground  state  are  forbidden.  Hence,  there  is  the  possibility  of  constructing 
atomic  frequency  standards  based  on  optical  photons  having  stabilities  character¬ 
istic  of  microwave  atomic  standards.  In  Section  2,  we  show  how  the  Raman 
process  results  in  a  narrow  resonance  whose  width  is  independent  of  the  inter¬ 
mediate  state  width. 


2.  THE  STIMULATED  RESONANCE  RAMAN  PROCESS 


Consider  the  energy  level  diagram  shown  in  Figure  1.  An  atom  in  level  1 
moving  at  a  speed  v  interacts  simultaneously  with  two  monochromatic  optical 
fields,  a  pump  field  Ep  at  frequency  Up  and  a  Stokes  field  Eg  at  frequency  ug. 
For  simplicity,  both  fields  are  assumed  to  be  zero  outside  the  interaction  region 
and  perfectly  uniform  in  the  interaction  region,  which  is  of  length  d  along  the 
atomic  velocity  vector,  such  that  the  transit  time  through  the  region  is  r  =  d/v. 
The  fields  at  the  position  of  the  atom  are  given  by 


-iu  t+i<4 

Ep(t)  =  (1/2)  ep  Eop  e  ?  P  +  e.  c. 


-iu  t+i$ 

Es(t>  -  (1/2)  eg  Eo|  e  8  s  +  c.c. 


(1) 


where  £p  and  eg  are  unit  polarization  vectors  for  the  waves  and  t  is  proportional 
to  the  distance  the  atom  has  penetrated  into  the  interaction  region.  In  order  that 
there  be  no  Doppler  shift,  the  E  field  wave  vectors  both  should  be  perpendicular 
to  the  atomic  velocity. 

The  dynamics  of  the  three-level  atom  is  described  by  Schrodinger's  equation 
with  a  phenomenological  decay  term  added 

»  If  -  H0  -  m  r*  .  (2) 

where  H  is  the  Hamiltonian  for  the  atom  and  Ib  the  sum  of  an  unperturbed  part  H° 
and  a  perturbation  V(t), 

H  =  H°  +  V(t)  ,  (3) 
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Figure  1.  Energy  Level  Diagram  Show¬ 
ing  Stimulated  Raman  Interaction  with 
Two  Optical  Fields 


and  r  is  a  phenomenological  relaxation  operator  describing  atomic  level  decay. 
Since  the  atom  is  much  smaller  than  an  optical  wavelength,  the  perturbation  can 
be  represented  by  a  dipole  interaction  term, 

V(t)  =  -  p  .  [Eft)  +  E  ft)]  ,  (4) 

where  n  is  the  atomic  electric  dipole  moment.  We  now  expand  the  wave  function 
0  in  unperturbed  atomic  states  ]m)  (m  =  1,  2,  3),  satisfying 

H°  |m)  =  Em  |m)  ;  (5) 

that  is, 


3 

0( t)  =  am(t)  exp  (-i  Emt/f>)  |m)  .  (6) 

m  =  l 

We  also  assume  that  mth  level  amplitude  decays  at  a  rate  7m/2,  so  that 

r  |m)  =  (ym/2)  |m)  .  (7) 

Substituting  Eqs.  (6)  and  (3)  into  Eq.  (2)  and  using  Eq.  (5),  Eq.  (7),  and  the 
orthogonality  of  the  eigenstates  |m),  (m|n)  =  a  we  obtain  the  equation  of 
motion  for  the  probability  amplitude  in  the  interaction  picture. 
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an  =  -(i/*>  2 

m  =  1 


iw„_,t 

e  Vnm(t)am-V2) 


(8) 


In  this  equation 


u>  =  (E 
nm  '  n 


(9) 


V  (t)  =  -u 
nm  t- nm 


[E  (t)  +  E  (t)] 
P  & 


(10) 


and  u  m  is  the  transition  moment  (off-diagonal  matrix  element  of  the  dipole  oper¬ 
ator);  the  superior  dot  denotes  time  differentiation.  The  square  modulus  of  a 

i  i  1 2  n 

gives  the  probability  of  being  in  the  state  |  n)  ,  =  |  |  .  and  the  phase  of  a 

gives  the  phase  shift  for  the  state,  whose  time  derivative  is  proportional  to  the 
energy  shift  of  the  level  due  to  the  interaction. 

To  obtain  Eq.  (8)  in  the  explicit  form  to  be  used,  we  substitute  from  Eqs.  (10) 
and  (1)  and  neglect  any  nonresonant  terms,  including  terms  varying  in  time  near 
twice  the  optical  frequency  such  as  exp  (±i(up  +  u21HI  or  exp  [±i(us  +  u23)t].  The 
result  is 

*  *  ‘4  +'Au  t 

&1  +  <7i/2)ai  =  (i/2h)  *»2i  EQp  e  p  P  a2  (11a) 


id  -iAw  t  i</>  -iAu  t 

a2  +  <?2/2>a2  =  (W  1*21  Eop  e  al  +  (i/2f,)  "23  Eos  e  a2 


*  *  -i*--i4wt 

a3  +  ^3/2^a3  =  M23  eos  e  a2 


(lib) 

(11c) 


where 


Awp  a  Up  -  w21 
^us  s  us  '  u23 

"21  =  ?p  ’  "21 

"23  =  es  '  " 23  • 


(12a) 

(12b) 
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We  assume  initially  at  time  t  =  t  when  the  atom  enters  the  interaction  region, 
there  is  no  population  in  states  |2)  and  |3), 


am(to)=5mlal(to'  • 


and  we  solve  the  equations  perturbatively  for  a^it),  to  lowest  possible  order  in  the 

dioole  moments  u  .  The  reader  who  is  not  interested  in  the  details  of  the  deriva- 
r  mn 

tion  may  jump  to  Eq.  (26). 

In  zero  order  only  a^  is  nonzero  and  satisfies  Eq.  (11a)  with  a2  =  0, 

®i(0>  +  (V2)ai<0>  =  0  •  (14) 

whose  solution  is 


Hence,  a^°\t)  decays  exponentially,  resulting  in  a  Lorentzian  line  broadening  of 
level  1  with  a  width  yr  In  a  similar  fashion.  y2  and  y3  will  result  in  widths  for 
levels  2  and  3.  respectively.  However,  as  discussed  in  Section  1,  we  will  ignore 
the  widths  of  levels  1  and  3,  that  is, 

7x<t  -  to)  ,  73(t  -  tQ)  «  1  ,  <16> 

and  assume  that  only  level  2  has  an  appreciable  width.  Hence,  Eq.  (15)  becomes 
simply 

ai<0\t)  a  ai(t0)  .  (17) 

To  first  order  in  p  ,  only  a.  is  nonzero  and  satisfies  Eq.  (lib)  with 
a  |  =  aj  '  and  a^  *  0, 


a2(l)  +  (72/2)a2(1>  =  (1/2*)  »2l  EQp  e^p  a/^t)  . 


Solving  Eq.  (18)  and  using  Eqs.  (16)  and  (17),  we  obtain  the  probability  amplitude 
for  state  2, 


a2(1>(t)  =  (1/2 tS)  M21  Eop  e'*P  -  iAUpl  -1 

f  -iAu_t  -y0 (t  — t  )/2 - 1 Aw  t  "I 

X  [e  P  -e  2  POJ  .  (19) 

If  we  look  for  time  differences  t  -  tQ  greater  than  the  reciprocal  linewidth  fj"1 
and  assume  that  the  frequency  up  is  within  a  linewidth  of  u>2^  —  that  is,  that  the 

Raman  process  is  resonant. 


y2{t  '  V  1  VA<Jp  ^  1 


(20) 


^  the  second  term  of  Eq.  (19)  can  be  neglected,  and  the  first  simplifies  to 


/i\  t) 

*2  <t)a(i*21Eop/^2)al(to)e  P  P 


(21) 


By  substituting  Eq.  (21)  for  a2  into  Eq.  (11c)  and  using  Eq.  (16),  we  obtain  an 
equation  for  a3,  second-order  in  p 


a3<2)  =  (i/2fi)  „23*E0*  e'l*s+lAWst  a2(1)(t)  , 


which  can  be  integrated  to  yield 

a3(2,(t)  =  [nR,  p  Vs  eiX»i(‘0V2i^2]  [e  ‘V  -  e  ‘V°]  .  (22) 

Here  fl0  and  nD  are  the  Rabl  frequencies  corresponding  to  the  pump  and 
n,  p  n,  s 

Stokes  fields, 

nR.  p  =  p21Eop/f' 

(23) 

nR,  s  =  m23  Eos/ft  • 

X  is  the  phase  difference  between  the  two  fields, 

X  ■  *p  -  *s  •  <24> 
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and 


AR  S  AuP  '  Aws  =  up  '  us  “  “31  •  <25) 

Defining 

g(x)  s  e  ix  sin  x/x  a  e  *x  sine  x  ,  (26) 

we  can  rewrite  Eq.  (22)  in  the  form*'® 


(2) 


<t)  =  -'%.pnR,s(t 


tQ)/2r2 


al(to}  e 


i(x-^Rt0) 


g{AR(t  -  to)/2j 


(27) 


The  transition  probability  to  level  3,  in  second-order  perturbation  theory,  is  then 
2  2 

W3<2)(t)  =  |a3(2)(T)|  =  knRpnR>sT/2T2)g(ART/2)|  .  (28) 

where  we  have  assumed  the  atom  enters  the  field  at  t  =0  and  remains  for  a  time 

o 

t.  From  Eq s.  (16)  and  (20),  we  can  summarize  the  conditions  for  the  validity  of 
Eq,  (28)  as 

yl  T  '  y3  T  «  1  «  ?2  r  ■  /  y2^Aup  '  <29) 

The  properties  of  the  stimulated  Raman  effect,  in  particular,  the  stimulated 

Raman  lineshape,  can  be  inferred  from  Eq,  (28).  It  is  clear  from  Eq.  (26)  that 
(2) 

W3  peaks  at  AR  =  0.  This  is  the  Raman  resonance  condition,  which  states 
explicitly 


“p'us  =  w31  •  (30) 

or  the  difference  between  the  pump  and  Stokes  frequencies  equals  the  atomic 

transition  frequency  between  initial  and  final  state.  Moreover,  as  the  optical 

frequency  difference  is  detuned,  the  atomic  response  falls  off  as  the  square  of 

o 

the  sine  function.  Since  sine  1.39  -  1/2,  it  follows  from  Eq.  (28)  that  the  Raman 
linewldth  (FWHM)  is 


rR  =  2(1. 39)/ xt  =  0.  885/t 


(31) 


10.  The  reader  who  has  not  followed  the  derivation  above  should  refer  to 

Eqs.  (23)  to  (26)  for  the  definitions  of  symbols  occurring  in  this  equation. 
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One  concludes  that  the  Raman  line  is  transit -time  broadened  and  insensitive  to 
the  broad  linewidth  of  level  2  (10  MHz  in  the  case  of  sodium).  In  a  very  real 
sense,  the  atom  responds  as  if  it  had  interacted  with  a  beam  of  microwaves  of 


frequency  u  =  w  -  u„. 

pps 


A  beam  of  sodium  atoms  emerging  from  an  oven  at  a 

4 

temperature  of  673K  (400 C)  has  most  probable  speed  v  =  9  X  10  cm/ sec. 

When  the  atoms  traverse  an  interaction  region  of  length  d  =  1  mm,  Eq.  (31)  pre¬ 
dicts  a  Raman  linewidth  of  Tp  =  800  kHz,  much  less  than  the  10  MHz  value  of 


By  increasing  the  length  of  the  interaction  region,  the  Raman  width  decreases 
in  inverse  proportion.  Nevertheless,  there  is  a  limit  to  the  amount  of  narrowing 
to  be  achieved  in  this  way,  if  for  no  other  reasons  than  that  a  finite  amount  of 
laser  power  is  available.  In  Section  3  we  examine  another  method  of  decreasing 
the  transit -time  linewidth  of  the  Raman  resonance. 


3.  THE  RAMSEY  METHOD  IN  STIMULATED  RAMAN  EMISSION 

The  experiment  to  detect  Raman  scattering  described  in  Section  2,  whose 

lineshape  and  width  are  given  in  Eqs.  (28)  and  (31),  is  analogous  to  Rabi's 

11  12  13  12 

resonance  method  in  magnetic  resonance.  '  In  1949  Ramsey  ’  proposed  a 
method  of  beating  transit  time  broadening  in  molecular  beam  magnetic  resonance 
experiments,  based  on  the  use  of  separated  oscillating  fields  or  separated  inter¬ 
action  regions.  The  Ramsey  method  has  also  been  applied  successfully  in  optical 
resonance  in  the  areas  of  saturated  absorption  spectroscopy  and  Doppler-free 
two-photon  absorption  spectroscopy,  starting  with  the  work  of  Baklanov  and 
others14,  15  in  1976. 

o 

To  apply  the  Ramsey  method  in  stimulated  Raman  scattering  spectroscopy, 
we  let  the  E  and  E  together  interact  with  the  atomic  beam  in  two  regions  of 
length  d  separated  by  a  much  larger  distance  L-d  (Figure  2).  The  distance  L 
corresponds  to  a  transit  time  T  =  L/v.  The  atomic  transition  to  level  3, 
described  by  Eq.  (27),  cam  occur  in  either  region  A  or  region  B,  the  upstream  and 


11.  Rabi,  1. 1.,  Zacharlas,  J.R.,  Millman,  S. ,  and  Kusch,  P.  (1938)  Phys.  Rev. 

5j3:318;  Kellogg,  J.M.  B.,  Rabi,  1. 1.,  Ramsey,  N.  F. ,  and  Zacharlas,  J.  R. 
TT939)  Phys.  Rev.  55:728. 

12.  See  also  Ramsey,  N.  F.  (1956)  Molecular  Beams  (Clarendon  Press.  Oxford). 

13.  Ramsey,  N.  F.  (1949)  Phys,  Rev.  76:996. 

14.  Baklanov,  Ye.V. ,  Dubetsky,  B.Ya.,  and  Chebotaev.  V.  P.  (1976)  A ppl.  Phys. 

9:171;  JU:201.  - *— 

15.  Bergqulst,  J.C.,  Lee,  S.H.,  and  Hall,  J.  L.  (1977)  Phys.  Rev,  Lett,  38:159. 
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Figure  2.  Experimental  Setup  for  Applying  Ramsey  Method  to 
Stimulated  Raman  Scattering  by  an  Atomic  Beam 


1 6 

downstream  regions  respectively.  In  perturbation  theory  the  transition  is  weak 
In  each  region,  and  the  transition  amplitude  small;  if  we  observe  the  atom  after  it 
has  traversed  both  regions,  the  total  transition  amplitude  is  the  coherent  sum  of 
the  individual  transition  amplitudes  for  each  region. 

Denote  the  amplitudes  for  undergoing  transitions  in  regions  A  and  B  by  sub¬ 
scripts  on  a^.  Let  t  =  t  =  0  when  the  atom  enters  region  A  and  a^(0)  =  1, 

Equation  <27)  then  yields 

a3.  A(2)(t)  -  -<«R.  pnR.  s  T/2f2>elXA  g(ARr/2)  .  (32) 

where  the  subscript  A  has  been  added  to  X  to  allow  the  phase  difference  of  the  two 

optical  waves  to  be  different  in  the  two  regions.  This  amplitude  then  propagates 

(2) 

freely  until  time  t  =  T  +  t  when  it  is  coherently  combined  with  a^  g  (T  +  t).  We 
assume  stronger,  yet  still  realistic,  conditions  on  and  y^  than  in  Eq.  (29), 

7lT.  73T  «  1  .  (33) 

so  that  ag  and  a^  do  not  change  when  they  propagate  freely.  Thus  we  have 

a3  a(2)(T  +  t)  =  a3  a(2)(t)  .  (34) 

Similarly,  to  find  the  amplitude  for  undergoing  a  transition  in  region  B,  we  let 
tQ  =  T  and  a^  g(T)  =  aj(0)  =  1,  obtaining 

16.  The  A  and  B  regions  defined  here  have  no  relation  to  the  A  (state  preparation) 
and  B  (detection)  regions  of  a  standard  atomic  clock.  In  fact,  the  regions 
called  A  and  B  here,  taken  together,  correspond  to  the  C,  or  interaction, 
region  of  an  atomic  clock. 
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g<ART/2) 


(35) 


a3,  B<2>(T  +  T>  *  ”<nR,  pnR%T/272>  ^ 


The  transition  probability  to  level  3  is  given  by 

2 

W3(2)(T  +  r)  =  |a3  a(2)(T  +  r)  +  a3  b(2)(T  +  t))  ;  (36) 

since  the  phase  factors  in  Eqs.  (32)  and  (35)  differ  in  a  frequency -dependent 

o 

manner,  we  expect  the  lineshape  to  differ  from  the  sine  x  lineshape  for  the  Rabi 
method  in  Eq.  (28).  Substituting  from  Eqs.  (32),  (34),  and  (35)  into  Eq.  (36),  we 
note  that  the  Raman  lineshape  in  the  Ramsey  method  is  given  by 


W3(2)(T  +  T) 


I^R,  p  nR,  s  T/?2^S^RT,/2^  W1  +  cos  (ART  +  *AB)1/2*  ’ 


(37) 


where  XAg  =  xA  "  Xg.  From  Eqs.  (31)  and  (29)  we  may  recapitulate  the  conditions 
under  which  this  expression  is  valid. 


yxT  ,  73T«  1  «  72t  ,  72/Aup  •  <38> 

The  Ramsey  lineshape  is  identical  to  the  Rabi  lineshape  except  for  the  factor 
in  curly  brackets.  As  the  field  frequency  is  tuned,  this  factor  varies  between  0 
and  1  at  a  rate  which  is  rapid  compared  to  the  factor  |  g|  2  because  T  »  t.  Hence, 
the  Rabi  lineshape  is  modulated  by  cosinusoidal  Ramsey  fringes  due  to  the  inter¬ 
ference  between  the  transition  amplitudes  in  the  two  regions.  Looked  at  in 
another  way,  the  first  (constant)  term  in  Eq.  (37)  represents  the  sum  of  the 
transition  probabilities  in  regions  A  and  B,  which  are  equal  in  this  approximation, 
while  the  second  (sinusoidally  varying)  term  is  a  cross  term  in  the  individual 
amplitudes. 

If  we  carefully  control  the  low-frequency  phase  difference  between  Ep  and  Es 
and  arrange  it  to  be  the  same  in  regions  A  and  B,  then  XAB  *  0.  In  this  case  the 
central  fringe  peaks  at  =  0,  and  the  width  (FWHM)  of  the  central  fringe  is 

Tr'  =  (2T)'1  .  (39) 

For  the  typical  sodium  atom  of  Section  2  and  L  =  20  cm,  this  results  In  rD'  *  2  kHz. 

n 

By  placing  the  interaction  regions  farther  apart,  we  reduce  the  Raman  linewldth 
in  the  Ramsey  method  yet  further.  Ultimately,  our  ability  to  reduce  the  linewldth 
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in  this  way  is  limited  only  by  practical  constraints  on  the  length  of  the  apparatus 
and  volume  of  the  vacuum  system;  by  the  mean-free  path  of  the  atoms  in  the  beam 
region;  and  by  the  drop  in  atomic  beam  flux  density  at  region  B  for  large  separa¬ 
tions,  combined  with  the  finite  amount  of  available  laser  power. 

Comparing  Eqs.  (37)  and  (31)  we  see  that  the  reduction  in  the  linewidth  by 
going  to  the  Ramsey  method  is  rR’/rR  ~  t/T  =  d/L.  For  the  typical  values  of 
d  and  L  assumed  above,  we  find  rR'/rR  =  1/200.  Since  the  interacting  atomB 
are  in  the  form  of  an  atomic  beam,  the  transition  probability  given  by  Eq.  (37)  or 
(28)  must  be  averaged  oVer  the  speed  distribution  in  the  atomic  beam.  However, 
this  averaging  will  not  alter  significantly  the  values  rRf  and  rR  of  the  linewidth, 
and  we  do  not  consider  its  effect  on  transit -time  broadening. 

In  this  section  and  Section  2  we  have  examined  the  Influence  of  one  broadening 
mechanism  —  transit -time  broadening,  a  result  of  the  short  time  which  the  atoms 
spend  in  the  laser  fields  —  on  the  Raman  process.  We  have  examined  transit -time 
broadening  in  the  Rabi  method  and  in  the  Ramsey  separated-field  method  where  its 
effect  is  diminished.  In  the  next  section  we  examine  other  sources  of  broadening. 


4.  BROADENING  MECHANISMS  IN  STIMULATED  RAMAN  EMISSION 

Among  the  other  possible  broadening  mechanisms  that  could  affect  the  stabil¬ 
ity  of  atomic  frequency  standards  based  on  stimulated  Raman  emission,  in  addition 
to  transit -time  broadening,  are  Doppler  broadening,  atomic  level  broadening, 

laser  jitter,  and  power  broadening  (or  optical  Stark  broadening).  Line  broadening 

2 

enters  into  the  frequency  stability  since  the  Allan  variance  a  (2,  t,  t),  the  usual 
measure  of  fractional  frequency  stability,  is  given  by 

o  =  (rT/vQ)/(S/N)  ,  (40) 

where  S/N  is  the  signal-to-noise  ratio,  the  transition  frequency,  and  rT  the 
total  linewidth.  Hence,  larger  llnewidths  require  larger  signal-to-noise  ratios 
to  obtain  a  given  fractional  frequency  stability.  We  consider  four  sources  of 
broadening  as  mentioned  in  the  foregoing  discussion. 

4.1  Doppler  Broadening 

When  a  Raman  transition  with  a  moving  atom  occurs,  the  resonance  frequency 
of  the  effective  microwave  field  which  the  atom  sees  is  shifted  from 
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to 


where  k  is  the  effective  microwave  wave  vector 


This  is  the  first-order  Doppler  effect.  By  keeping  k  orthogonal  to  v  we  can  make 

— 

the  Doppler  frequency  shift  equal  to  2ero.  However,  with  an  atomic  beam  it  is 
impossible  to  do  this  for  all  atoms,  since  the  beam  has  a  nonzero  divergence. 

There  is  a  residual  shift  due  to  the  radial  component  of  the  beam  velocity,  which, 
for  small  beam  divergence  6  takes  on  a  maximum  value 

v  =  u0 
r  o 

When  the  velocity  average  is  performed,  the  net  effect  of  all  the  shifts  is  to 
broaden  the  Raman  transition,  leading  to  a  Doppler  linewidth 

rD  =  [(ln2)1/2/r]  k(juflQ  .  (42) 

For  the  1.  77  GHz  hyperfine  transition  of  sodium  and  an  atomic  beam  divergence 
of  1  mrad,  one  easily  finds  rD  =  5  Hz.  Hence,  this  source  of  broadening  is 
quite  negligible  compared  to  transit -time  broadening. 

4.2  Atomic  Level  Broadening 

A  second  potential  source  of  broadening  is  the  linewidths  of  the  initial  and 
final  states  1  and  3  of  the  atom.  These  were  denoted  and  y3  in  the  Eqs.  (11)  and 
quite  properly  neglected  in  Eq.  (16).  The  radiative  lifetimes  of  levels  1  and  3  are 
extremely  long,  and  there  is  no  collision  broadening  of  the  levels  in  an  atomic 
beam. 

4.3  Later  Jitter 

The  Raman  resonance  can  be  broadened  not  only  by  the  widths  of  atomic  states 
l  and  3,  but  also  by  the  spectral  width  of  the  light  incident  on  the  atom.  If  we 
remove  the  requirement  that  the  laser  light  fields  be  monochromatic,  then 
Eqs.  (1)  are  replaced  by 
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-iUpt-i  J  6u)p(t'>dt,+i^p 

Ep(t)  =  <l/2)epEope  °  +  c.  c. 

(43) 

t 

-iust-i  5  6us(t')dt'+i(J>g 

Eg(t)  =  (l/2)esEQse  0  +c.c. 


Here  the  instantaneous  frequencies  6up  and  6ug  are  random  functions  of  time 

accounting  for  frequency  jitter  of  the  two  light  beams.  The  difference  of  the  phase 

of  E_  and  Eg  in  Eq.  (43)  enters  the  calculation  of  the  Raman  transition  amplitude 

a^l2*.  Since  6wp  and  6ug  are  statistically  independent  processes,  the  variance  of 

the  difference  frequency  is  the  sum  of  the  variances  of  6 u  and  6u  and  the  rms 

p  s 

contribution  to  the  linewidth  is 


rT  -  (<«0  +  <«</>) 


1/2 


(44) 


2  I/2 

where  the  angular  brackets  denote  statistical  averaging.  Since  (6cjp)  can  be 

hundreds  of  kilohertz,  laser  jitter  would  result  in  an  unacceptably  large  linewidth 

completely  dominating  the  transit -time  width. 

Nevertheless,  there  is  a  rather  simple  way  to  eliminate  the  effect  of  laser 

jitter,  and  that  is  to  arrange  for  the  jitter  in  the  two  laser  fields  to  be  correlated. 

This  is  accomplished  by  locking  the  two  fields  to  a  common  reference  cavity  or, 

more  simply,  if  the  frequency  difference  is  in  an  accessible  region  such  as  the 

microwave  region,  by  deriving  w  from  the  laser  frequency  u  through  electro- 

1 7P  s 

optic  or  acousto-optic  shifting.  The  acousto-optic  shifting  in  a  Bragg  diffraction 
cell  is  particularly  efficient  and  is  convenient  to  carry  out  at  the  sodium  transition 
frequency  of  1.  77  GHz.  Letting  (t)  be  the  instantaneous  laser  frequency  due  to 
jitter,  we  find 


ug(t)  =  Uj  +  iWj(t) 
u  <t)  -  Uf  +  Wp  +  (t)  . 


17.  Lelby,  Jr.,  C.C.,  and  Ezekiel,  S.,  Laser  Excited  Molecular  Beam  Time  and 
Frequency  Standard,  allowed  U.S.  patent  application  No.  134,358,  filed 
27  March  1980. 


21 


When  the  difference  frequency  ug  -  up  is  formed,  the  laser  jitter  drops  out  of  the 
problem, 

ws(t)  ‘  wp(t)  =  %i  * 
independent  of  6uj(t). 

To  determine  the  effect  of  laser  jitter  in  a  more  precise  manner,  we  can  re¬ 
derive  equations  of  motion  for  the  probability  amplitudes  an(t)  as  we  did  in 
Section  2,  but  using  Eqs.  (43)  rather  than  Eqs.  (1)  for  the  optical  fields.  The 
result  is  Eqs.  (11)  with  replaced  by  ^  +  6^.(t)  (i  =  p,  s),  where  we  have  defined 

t 

6*i<t)  =  -  Jfiujit')  dt*  .  (45) 

o 

If  we  carry  out  a  derivation  analogous  to  that  of  Eq.  (27),  we  find  that  Eq.  (19)  is 
replaced  by 


( 1 )  ^  # 

a2*  (t)  =  (i/2)nH,p  e  P  al(to)  /  exp  *  t'V2  *  i£*>pv  +  i**p<t')]  dt* 


(46) 


while  Eq.  (22)  is  replaced  by 


-lADt' 


a3(2>(t)  =  -U/4)nR>pn^8ai(to)e1X/dt.e'  « 


t' 

X  f  dt"  exp  {-T2(f  -  t")/2  +  i[5*p(t")  -  6#s(t')l  }  .  (47) 

*o 

In  deriving  Eq.  (47)  we  have  used  the  approximations  in  Eqs.  (16)  and  (20). 

Assume  now  that  the  phases  of  the  pump  and  Stokes  waves  are  correlated,  so 

that 


22 


* 


«*8<t)  =  6*p(t)  ■  6*<t) 

and  (48) 

4ws(t)  =  «up(t)  *  6ut  (t)  . 

Then,  for  very  slowly  varying  phases,  we  can  expand  6 <fi  in  a  Taylor  series, 

6*<t")  -  6^(t')  et  -  t')  =  -6ut(f)(t"  -  t> )  ,  (49) 

and  integrate  Eq.  (47)  to  obtain 


a3<2)(t)  =  -  <nR  pnR*aaj<t0)  eiX/2y2)  f  dt»  e  [1  -  2i6U|  <t')/y2]  _1  . 

t 

o 

(50) 

It  follows  from  Eq.  (50)  that  if  6u(  «  y2>  laSer  jitter  has  a  negligible  effect  on  the 
Raman  transition. 

The  final  integration  in  Eq.  (50)  cannot  be  carried  out.  Moreover,  the  series 
expansion  of  Eq.  (49)  is  valid  only  in  the  limit  Y2TC  »  1,  where  rc  is  the  correla¬ 
tion  time  of  the  frequency  fluctuations.  An  approach  which  is  not  restricted  to  this 
limit  consists  in  averaging  the  random  process  £u((t)  over  an  appropriate  proba¬ 
bility  density  functional.  To  be  correct,  we  should  average  the  transition 
probability18  Wg(2)  =  |a3(2)i2;  however,  we  obtain  qualitatively  correct  results  by 
averaging  the  transition  amplitude  in  Eq.  (47), 


(2) 


t  t> 


*  iv  f  t  -Y2(t,-t")/2 

=  -<l/4)nR >pfiR  gai(to)  elX  J  dt'  e  J  dt"  e  Q[l;t",t'], 


where 


(51) 


Q(X(t);t",  t'] 


dt  X(t)  6ut  (t) 


(52) 


18.  Picard,  R.H.,  and  Willis,  C.  R.  (to  be  published). 
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is  the  characteristic  or  moment-generating  functional  for  the  random  process 
6u(  (t). 

For  typical  sources  of  laser  jitter,  such  as  vibration  of  a  laser-cavity  mirror, 
it  is  reasonable  to  assume  that  6uf(t)  is  a  stationary,  zero-mean,  Gaussian  proc¬ 
ess;  then  Q  is  determined  solely  by  its  second  moment  or  correlation  function 


C(t1  -  t2)  a  <6«j (tj)  6ut( t2)> 
19 

through  the  relation 


(53) 


Ql\(t);t",  t'l  =  exp  /-(1/2)  f  dtj  j  dt2  C(t1  -  tg)  X(tj)  X(t2)>  . 

I  t"  t"  J 


(54) 


After  setting  X(t)  =  1,  using  statlonarity  in  the  form  C(-r)  =  C(r),  and  changing  the 
order  of  integration,  one  notes  that  Eq.  (54)  yields 


Qlljt",  t']  =  exp 


dr(t' 


t)C(t)^  , 


(55) 


Notice  that,  with  X(t)  =  l.Q  is  dependent  only  on  the  difference  of  its  time  argu¬ 
ments  t'  -  t",  so  that  we  may  represent  it  by  Q(l;t'  -  t"] . 

Equations  (51)  and  (55)  comprise  the  complete  solution  of  our  problem.  Oice 

the  correlation  function  C(t)  is  specified,  we  can  compute  the  characteristic 

(2)  , 

functional  Q  and  use  it  to  determine  the  transition  amplitude  a^  ,  To  determine 
the  effect  of  jitter  on  the  width  of  the  atomic  response,  we  can  change  variables 
of  integration  in  Eq.  (51)  to  t<  a  t"  -  tQ  and  r"  *  t'  -  t".  then  Interchange  the  order 
of  integration,  with  the  result 


t-t 


j<2)  =  -(1/4)  f  dr" 


-I(y2/2)+lARli 


t-t  -r" 


QIIjt"] 


/ 


dT'  e 


-iART’ 


(56) 

Since  r"  is  limited  to  lt  can  be  l8norert  compared  to  t  -  tQ,  by  Eq.  (20), 

and  the  inner  integral  has  the  value  (t  -  tQ)g(AR(t  -  tQ)/2] .  If  we  restrict  ourselves 


19.  Middleton,  D.  (1960)  Statistical  Communication  Theory.  McGraw-Hill, 
New  York,  pp.  372^37 T. 
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to  detunings  satisfying  AR  «  y2.  a  very  weak  restriction,  we  may  also  ignore 
in  the  exponent  of  the  t"  integral,  so  that  finally  a3  is  given  by  Eq.  (27)  with  the 
right  member  multiplied  by  the  factor 


-y2r"/2 


(57) 


Further  details  on  the  influence  of  laser  jitter,  including  the  important  limit- 

18 

ing  cases  of  slow  modulation  and  fast  modulation,  will  be  presented  elsewhere. 
However,  important  qualitative  conclusions  on  laser  jitter  of  arbitrary  strength 
and  time  scale  can  be  drawn  from  Eq,  (57).  Since  the  factor  in  Eq.  (57)  is  gen¬ 
erally  smaller  than  unity,  jitter  will  lead  to  a  small  reduction  of  the  transition 
amplitude.  However,  because  Eq.  (57)  is  real  and  independent  of  AR,  jitter  will 
not  affect  the  atomic  lineshape  or  linewidth  at  all.  Hence,  laser  frequency  fluctua¬ 
tions,  like  the  broad  width  of  level  2,  do  not  affect  the  Raman  linewidth,  as  long 
as  we  correlate  the  jitter  at  the  two  frequencies  in  some  fashion. 

4.4  Power  Broadening 

Finally,  we  mention  the  influence  of  power  broadening.  When  the  field  Ep  is 

applied  to  the  atom  near  resonance  between  levels  1  and  2,  both  levels  split  by  an 

amount  equal  to  the  Rabi  frequency  flR  p  (optical  Stark  effect).  If  the  splitting  is 

small  compared  to  the  level  widths,  it  is  not  resolved  and  appears  as  a  broadening 

2 

of  the  level  by  an  amount  nR  p  /2ry^,  when  y Aup  «  y 2.  This  process  is  illus¬ 
trated  in  Figure  3(a).  The  interpretation  of  this  broadening  is  that  level  1  has 
acquired  a  width  due  to  the  coherent  superposition  of  state  1  and  the  short-lived 
state  2  induced  by  the  optical  field.  In  a  similar  fashion,  the  Stokes  field  broadens 

levels  3  and  2,  with  the  result  that  the  Raman  transition  from  1  to  3  in  the  Rabi 

20 

method  will  have  a  width  due  to  power  broadening  [Figure  3(b)]  given  by 

rp  =  (\p2  +  f!R,>^2  *  (58) 

Here  it  has  been  assumed  that 


Ti.Tg.^p,  Aus  «T2  •  (59) 

•18 

For  dipole  moments  typical  of  sodium,  p ^  an<*  P32  s  ®  Debye  (1  Debye  =  10  esu), 
and  light  intensities  I  =  Ig  =  50  pW/cm^,  we  find  Tp  =  250  kHz.  This  is 

20.  See  Bloembergen,  N.,  and  Levenson,  M.D.  (1976)  in  High-Reaolution  Laser 
Spectroscopy.  K.  Shimoda,  Ed.,  Springer-Verlag,  Berlin,  p.  315. 
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comparable  to  the  transit -time  broadening  rR  and  shows  that  one  must  limit  one¬ 
self  to  weak  laser  fields  to  avoid  making  power  broadening  the  dominant  broadening 
mechanism. 


(a)  PUMP  FIELD  ALONE  (•>)  PUMP  and  STOKES  FIELD 

Figure  3.  Effects  of  Optical  Stark  Splitting  on  Atomic  Level  Widths:  (a)  with  pump 
field  alone  present:  (b)  with  both  pump  and  Stokes  fields  present 


In  the  Ramsey  method  there  is  a  dramatic  reduction  in  power  broadening 
from  that  calculated  above  for  the  Rabi  method.  To  see  this  effect  qualitatively, 
we  recall  from  Eq.  (37)  that  the  Ramsey-method  lineshape  is  a  product  of  cosinu¬ 
soidal  interference  fringes  and  an  envelope  given  by  the  Rabi  lineshape.  The 
envelope  is  power  broadened  as  in  the  previous  paragraph.  The  fringes,  which 
arise  due  to  the  transit  of  the  atoms  through  the  dark  region  between  A  and  B.  are 
not  broadened  directly.  However,  when  the  indirect  effect  of  the  envelope's 
broadening  on  the  width  of  the  central  fringe  is  calculated,  it  is  found  that  the 
latter  is  broadened  by  a  factor  of  order  (d/L)^  times  the  value  in  Eq.  (58).  For 
the  parameters  assumed  in  Section  3,  this  reduction  factor  is  of  order  2  X  10  , 

and  the  resulting  width  is  of  order  10  Hz. 

In  passing,  we  note  that,  not  only  the  Unewidth,  but  also  the  slgnal-to-noise 
ratio,  is  laser  power  dependent.  The  latter  goes  up  with  the  square  root  of  the 
laser  power,  if  the  detection  proce  it-  photon  shot  -noise-limited.  However,  this 
increase  is  not  sufficient  to  overcome  the  increase  in  linewidth.  The  over-all 
result  is  that  the  Allan  variance  goes  up  with  power,  as  is  obvious  from  Eq.  (40). 


4.S  Implication*  for  Frequency  Standard* 

It  is  therefore  to  be  concluded  that,  with  proper  precautions,  none  of  the 
aforementioned  broadening  mechanisms  compete  seriously  with  transit -time 
broadening  as  an  important  constituent  of  the  Raman  linewidth  in  the  Ramsey 
method.  In  particular,  we  have  shown  how  the  effect  of  laser  jitter  may  be  com¬ 
pletely  eliminated  by  correlating  the  frequency  jitters  of  the  pump  and  Stokes  beams. 

5.  LINES  HI  FTS  IN  STIMULATED  RAMAN  EMISSION 

In  addition  to  the  broadening  introduced  by  the  various  physical  phenomena 
discussed  in  Section  4,  most  of  these  phenomena  also  introduce  a  shift  in  the 
position  of  the  Raman  line.  Determination  of  these  shifts  is  especially  impor¬ 
tant,  since  they  introduce  systematic  errors  into  the  frequency  standard.  Unlike 
line  broadening,  no  amount  of  Increase  in  slgnal-to-noise  ratio  will  decrease  the 
effect  of  the  lineshift  on  the  accuracy  of  the  standard. 

The  main  sources  of  Raman  lineshifts  are  first-order  and  second-order 
(transverse)  Doppler  effect,  differences  in  microwave  phase  <t>u  3  4>p  ~  4S  between 
the  two  interaction  regions  in  the  Ramsey  method,  and  power-dependent  atomic 

level  shifts  (light  shifts).  The  subject  of  lineshifts  is  important  enough  to  merit 

18 

attention  on  its  own.  A  separate  paper  will  be  published. 
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